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ABSTRACT 

We use a sample of 151 local non-blazar AGN selected from the INTEGRAL all¬ 
sky hard X-ray survey to investigate if the observed declining trend of the fraction of 
obscured (i.e. showing X-ray absorption) AGN with increasing luminosity is mostly 
an intrinsic or selection effect. Using a torus-obscuration model, we demonstrate that 
in addition to negative bias, due to absorption in the torus, in finding obscured AGN 
in hard X-ray flux limited surveys, there is also positive bias in finding unobscured 
AGN, due to Gompton reflection in the torus. These biases can be even stronger tak¬ 
ing into account plausible intrinsic collimation of hard X-ray emission along the axis 
of the obscuring torus. Given the AGN luminosity function, which steepens at high 
luminosities, these observational biases lead to a decreasing observed fraction of ob¬ 
scured AGN with increasing luminosity even if this fraction has no intrinsic luminosity 
dependence. We find that if the central hard X-ray source in AGN is isotropic, the 
intrinsic (i.e. corrected for biases) obscured AGN fraction still shows a declining trend 
with luminosity, although the intrinsic obscured fraction is significantly larger than 
the observed one: the actual fraction is larger than ~ 85% at erg s“^ (17- 

60 keV), and decreases to ^60% at L ^ 10'^'^ erg s“^. In terms of the half-opening 
angle 0 of an obscuring torus, this implies that 0^30° in lower-luminosity AGN, and 
0 ^ 45° in higher-luminosity ones. If, however, the emission from the central SMBH is 
collimated as dL/dfl oc cos a, the intrinsic dependence of the obscured AGN fraction 
is consistent with a luminosity-independent torus half-opening angle 0 ~ 30°. 

Key words: galaxies: active - galaxies: nuclei - galaxies: Seyfert. 


1 INTRODUCTION 


A lot of recent studies based on X-ray and hard X- 
ray extragalactic surveys have demonstrated that the 
fraction of X-ray absorbed (hereafter referred to as 
obscured) active galactic nuclei (AGN) decreases with 
increasing observed X-ray luminosity, at least at ^ 
10"^^ erg s ~^, both in the local ~ 0) and hi g h-redshift 
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lier evidence. iLawrence fc ElvislIlQsj ). This might indicate 
that the opening angle of the (presumably) toroidal obscur¬ 
ing structure - the key element of AGN unification schemes 
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- increases with AGN luminosity, for example due to feed¬ 
back of the central supermassive black hole (SMBH) on the 
accretion flow. 

Could the observed luminosity dependence of the 
obscured AGN fraction arise due to selection effects? 
This question has bee n occasionally raised before (e.g. 
iMavo &: Lawrencg|2013l) and is prompted by the fact that 
even hard X-ray (> 10 keV) surveys, which are usually 
flux (or signal-to-noise ratio) limited, should be biased 
against detection of Compton-thick AGN, i.e. objects viewed 
through absorption column density Ah > 10^^ cm“^, let 
alone X-ray surveys at energies below 10 keV which must 
be biased against even Compton-thin obscured sources. 
Due to this detection bias, the observed fraction of ob¬ 
scured AGN is expected to be lower than the intrinsic 
fraction of such objects. Furthermore, this effect may de¬ 
pend on luminosity, somehow reflecting the shape of the 
AGN luminosity f uncti on (LF). In fact, as discussed by 
[Lawrence fc Elvi3 (l2010l ). some mid-infrared selected, radio 
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selected and volume-limited AGN samples do not demon¬ 
strate any clear luminosity dependence of the proportion of 
type 1 (i.e. containing broad emission lines in the optical 
spectrum) and type 2 AGN. 

Although there have been previous attempts 
jUedae^J l2003l : iMalizia et al.l l2009l : iBurlon et al.l 1201 ll : 
lUeda et alj |2014| ~) to take into account detection biases 
when estimating the space density of obscured AGN based 
on hard X-ray surveys, they were, in our view, not fully 
self-consistent and/or used too small samples of hard X-ray 
selected AGN. It is our goal here to improve on both of 
these aspects. 

The purpose of the present study is to i) evaluate the 
impact on the observed hard X-ray LF and observed lumi¬ 
nosity dependence of the obscured AGN fraction of the neg¬ 
ative bias for obscured AGN discussed above and a positive 
bias that we demonstrate likely exists for unobscured AGN, 
and ii) reconstruct the intrinsic dependence of the fraction 
of obscured AGN on luminosity in the local Universe. Our 
treatment is based on a realistic torus-like obscuration model 
and makes use of the INTEGRAL/IBIS 7-year (2002-2009) 
hard X-ray survey of the extragalactic sky. Our sample con¬ 
sists of ~ 150 local (2;^0.2) Seyfert galaxies and is highly 
complete and reliable. Although there are now significantly 
larger hard X-ray selected samples of local AGN, based on 
additional observations by INTEGRAL/IBIS and especially 
by Swift/BAT, they currently suffer from significant incom¬ 
pleteness as concerns identification and absorption column 
density information. Most importantly, our sample is large 
enough to contain a significant number, 17, of heavily ob¬ 
scured (Nh ^ 10^"^ cm~^) AGN, for which we use as much 
as possible Nh estimates based on high-quality hard X-ray 
spectral data, in particular from the NuSTAR observatory, 
which has recently been systematically observing AGN dis¬ 
covered in the Swift/BAT and INTEGRAL/\B\S hard X-ray 
surveys. 

We adopt a AC DM cosmological model with flm = 0.3 
and Ho = 70 km s“^ Mpc“®. 


2 THE INTEGRAL AGN SAMPLE 


We use the catalogue of sources dKrivonos et al. I l2010bl) 
from the INTEGRAL/\B\S 7-year all-sky hard X-ray sur¬ 
vey (h ereafer, the INTEGRAL 7-year survey. iKrivonos et al.l 
l2010al '). To minimise possible biases in our study of the local 
AGN population due to remaining unidentified INTEGRAL 
sources and objects with missing distance and/or X-ray ab¬ 
sorption information, we exclude from the consideration the 
Galactic plane region (|6| < 5°). The catalogue is composed 
of sources detected on the time-averaged (December 2002 
- July 2009) 17-60 keV map of the sky and is significance 
limited (So). The corresponding flux limit varies over the 
sky: /det < 2.6 x 10“^^ (< 7 x 10“^^ erg s“^ cm“^) for 50% 
(90%) of the extragalactic (|6| > 5°) sky (see Fig. |T|. 

The main properties of the INTEGRAL 7-year sur¬ 
vey and of the correspond i ng cata logue of sources were de- 
scribed by Krivonos et all (I2010al f5l. Using this catalogue, 
ISazonov et al.l ( 2010l l rnade preliminary estimates of the 
hard X-ray LF of local AGN and the dependence of the 
obscured AGN fraction on luminosity. Subsequent follow-up 
efforts by different teams have resulted in additional identi- 



Limitlng flux (17 — 60 keV), erg s"' cm“^ 

Figure 1. Cumulative fraction of the extragalatic (|6| > 5°) sky 
as a function of flux limit in the INTEGRAL 7-year survey. 


fications, classifications, distance measurements and X-ray 
absorption column estimates for many INTEGRAL sources, 
which has significantly improved the quality of the cata¬ 
logue, as detailed below. 

The final sample used here consists of 151 non-blazar 
(i.e. Seyfert-like) AGN (see Table [AT] in Appendix M. with 
blazars (15 in total) being excluded from the analysis. The 
sample is highly complete, as there are only 4 sources at 
|6| > 5° from the INTEGRAL 7-year catalogue that remain 
unidentified. Moreover, all of our AGN have known distances 
and reliable estimates of their absorption columns based on 
X-ray spectroscopy. As illustrated in Fig. [2j our sample is 
mostly local, with 146 out of the 151 objects being located 
at 2 < 0.2, and spans about 5 decades in (observed) lumi¬ 
nosity, from Lobs ~ 10'^'^ to ~ 10“^® erg (hereafter, all 
luminosities are in the 17-60 keV energy band, unless spec¬ 
ified otherwise). 

We note that although we used the most up-to-date 
information from the NASA/IPAC Extragalactic Database 
(NED) and recent literature to remove blazars from our 
AGN sample, we cannot rule out that some of our objects 
have blazar-like properties, i.e. their observed hard X-ray 
emission contains a significant contribution from a relativis¬ 
tic jet. The most suspicious in this respect are objects clas¬ 
sified as broad-line (i.e. presumably oriented towards us) ra¬ 
dio galaxies. There are 6 such AGN in our sample: 3C 111, 
3C 120, Pic A, 3C 390.3, 4C -674.26 and S5 2116-681. All 
of them have Lobs > 10'*'* erg s“* (but < 10*® erg s“*), i.e. 
belong to the high-luminosity part of the sample. However, 
the total number of objects with Lobs > 10** erg s“* is much 
larger: 42. This suggests that possible incomplete filtering of 
the sample from blazars is unlikely to significantly affect the 
results and conclusions of this work. 


2.1 Absorption columns, heavily obscured AGN 

For the purposes of this study it is important to have 
maximally complete and reliable information on the X- 


© 0000 RAS, MNRAS 000, 000-000 



























Does obscured AGN fraction depend on luminosity? 3 


Table 1. Heavily obscured AGN from the INTEGRAL 7-year survey 


Object 

D 

Mpc 

-^obs 

erg 

Ah 

cm“2 

Reference for A^h 

SWIFT J0025.8-1-6818 

52.0 

3.2 X 

1042 

> 1023 

NuSTAE. fKrivonos et al. 2ni5j 

NGC 1068 

12.3 

3.8 X 

1041 

> 1023 

NuSTAR fBauer et al. 2014') 

NGC 1194 

59.0 

6.6 X 

1042 

~ 1024? 

XMM-Newton (below 10 keV. Greenhill. Tilak. X. Madejski 2008") 

CGCG 420-015 

129.2 

3.3 X 

1043 

> 1025 

NuSTAR (Krivonos et al. 2015^ 

MRK 3 

58.6 

2.8 X 

1043 

1024 

Suzaku (Ikeda^ Awaki. Kr. Terashima 2009j 

IGR J09253-I-6929 

172.6 

4.5 X 

1043 

> 1024? 

low X-ray/hard X-ray flux ratio (Swift/XRT+INTEGRAL/IBIS) 

NGC 3081 

28.6 

4.4 X 

1042 

1024 

Suzaku (^guch^^l^^OH') 

NGC 3281 

46.3 

1.2 X 

1043 

2 X 1024 

BewoSAX Q^ignal^^ Comastri 2002") 

ESO 506-G027 

109.5 

5.8 X 

1043 

1024 

Suzaku (Winter et al. 2009b'l 

NGC 4939 

34.7 

2.3 X 

1042 

> 1025? 

BevvoSAX (Maiolino et al. 1998). but varied to N-a = 1.5 X 1023 cm~2 
(XMM-Newton, below 10 keV. Guainazzi et al. 2005aj 

NGC 4945 

3.4 

2.6 X 

1041 

4 X 1024 

NuSTAR (Puccetti et al. 2014: Brierhtman et al. 20151. Suzaku (Yaaoob 2012') 

IGR J14175-4641 

348.3 

1.6 X 

1044 

> 1024? 

low X-ray/hard X-ray flux ratio (Swift/XRT+INTEGRAL/lBlS) 

NGC 5643 

11.8 

1.7 X 

1041 

> 1023 

NuSTAR (Krivonos et al. 20151 

NGC 5728 

24.8 

3.2 X 

1042 

2 X 1024 

NuSTAR. (Krivonos et al. 20151 

IGR J14561-3738 

107.7 

1.6 X 

1043 

~ 1024 

Chandra-\-INTEGRAL/IBlS (Sazonov et al. 20081 

ESO 137-G034 

33.0 

1.8 X 

1042 

> 1025 

Suzaku (Gomastri et al. 20101 

NGC 6240 

107.3 

5.8 X 

1043 

2.5 X 1024 

NuSTAR (Krivonos et al. 20151 



Figure 2. Observed hard X-ray (17-60 keV) luminosity vs. red- 
shift for non-blazar AGN from the INTEGRAL 7-year survey. 
Filled circles, empty squares and stars denote unobscured, lightly 
obscured and heavily obscured objects, respectively. 


ray absorption columns, A^h, of the studied AGN. Our 
starting source of such information is our previous pa¬ 
pers on the INTEGRAL/ISIS survey dSazonov et al.l[2007l . 
I2OI2I) as well as on the RXTE (3-20 keV) slew survey 
i Sazonov fc Revnivtsevir2004l ') . but we have updated the Ah 
estimates in all cases where it was necessary and possible 
(see Table IXTl) . 

For unobscured and lightly obscured (Ah < 10^'* cm“^) 
sources. X-ray spectroscopy at energies below 10 keV is usu¬ 
ally sufficient for evaluating Ah. Such data do exist for all 
of our sources and in most cases there are reliable published 
Ah values, which we adopt. Furthermore, if the absorption 


column is less than 10^^ cm“^, we adopt Ah = 0 and con¬ 
sider such sources unobscured. 

However, absorption column estimates based on X-ray 
data below 10 keV become unreliable for strongly absorbed 
sources, having Ah ^ 10^"^ cm~^. In such cases, we prefer 
to use results from hard X-ray (above 10 keV) spectroscopy, 
whenever possible. Specihcally, our preference list of instru¬ 
ments is headed by NuSTAR - the unique focusing hard 
X-ray telescope, followed by Suzaku and then by all other 
currently operating or previously flown hard X-ray missions. 

For hve of the heavily obscured (Ah > lO^'^ cm ob¬ 
jects and candidates, we carrie d out our own analysi s of pub¬ 
licly available NuSTAR data (|Krivonos et aklboish . Specif¬ 
ically, we htted the spectra by a sum of a strongly absorbed 
power-law component (with a high-energy cutoff) and a 
disk-reflection continuum modelled with pexrav in XSPEC. 
The NuSTAR spectra of SWIFT J0025.8-h6818, CGCG 420- 
015 and NGC 5643 are consistent with being fully reflec¬ 
tion dominated (i.e. dominated by Gompton-scattered con¬ 
tinuum), and so we prescribed Ah > 10^® cm~^ to them. 
The other two objects, NGG 5728 and NGG 6240, along with 
strong reflection demonstrate a significant contribution from 
the primary component suppressed by intrinsic absorption 
at the level of Ah ~ 2-2.5 x 10^'* cm“^. More physically 
motivat ed AGN torus models confirmed this qualitative re¬ 
sult (see iKrivonos et al.ll2015l for details). Our derived spec¬ 
tral parameters for GGCG 420-015, NGC 5643, NGC 5728 
and NGC 6240 ar e consistent with pie-NuSTAR estimates 

20nl. iMatt eralll20ll 


for these objects (lSeyergnini efLaI 


IComastri et al.ll2010l . IVignati et al 


199£ 


, respectively). 


In total, our sample consists of 67 unobscured (Ah < 
10^^ cm“^) and 84 obscured (Ah ^ 10^^ cm~^) AGN, in¬ 
cluding 17 heavily obscured (Ah 5: lO^'^ cm“^) ones. 

Table [T] provides key information about our heav¬ 
ily obscured AGN. For 7 of these, there are reliable 
Ah estimates or evidence that the source’s spectrum is 
reflection-dominated (in which case we adopt that Ah > 
10^® cm“^) from NuSTAR observations. All but one 
(IGR J14561—3738) of the remaining 10 objects are either 
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planned to be observed by NuSTAR soon or have already 
been observed by this telescope but the data are proprietary 
at the time of writing. However, for most of these sources 
there exists fairly reliable information from other hard X-ray 
missions indicating that A^h ^ 10^“* cm“^ - see Table [T] 

Three of the objects included in our sample of heav¬ 
ily obscured AGN are currently candidates rather than 
firmly established representatives of this class: the quoted 
value Ah 10^"' cm~^ for NGC 1194 comes from X- 
ray data below 10 keV, whereas the presence of Ah > 
10^^ cm“^ absorption columns in IGR J09253-I-6929 and 
IGR J14175—4641 is strongly suggested by very low (;^0.01) 
X-ray/hard X-ray flux ratios that we find for them from 
S’uij/t/XRT and INTEGRAL/\YAS> data. Note that we ini¬ 
tially used the same argument to regard another source 
from this sample, SWIFT J0025.8+6818, as a likely heav¬ 
ily obscured AGN, and it indeed proved to be such once we 
analysed NuSTAR data. In the analysis below, we assume 
that Ah = 3 x 10^^ cm~^ for both IGR J09253+6929 and 
IGR J14175-4641. 

The most difficult case is that of NGC 4939, 
which manifested itself as a reflection-dominated source 
(Ah > 10^^ cm~^) d uring BeppoSAX observations in 1997 
IIMaiolino et al.l[l998l '). but was found to be in a Compton- 
thin state , with Ah ^ 1.5 x 10^® cm“^, by XMM-Newton 
in 2001 IIGuainazzi et al.l l2005al ). We nevertheless treat 
NGC 4939 as a reflection-dominated source in our analy¬ 
sis, in part because the hard X-ray flux measured by INTE¬ 
GRAL for this source is similar to that measured by Bep¬ 
poSAX but lower than the flux inferred from the XMM- 
Newton observation and so INTEGRAL may have caught 
the source in a state similar to that revealed by Bep¬ 
poSAX. Generally, we adopt Ah = 10^® cm“^ for reflection- 
dominated sources (there are in total 6 such objects) in our 
analysis, although in reality the column density in such ob¬ 
jects may be even higher, say Ah ~ 10 ^® cm~^. 

We have thus obtained a fairly large and high-quality 
(in terms of information on intrinsic obscuration) sample of 
heavily obscured AGN. The high completeness and reliabil¬ 
ity of this sample are crucial for our analysis below. 


3 OBSERVED PROPERTIES OF LOCAL AGN 

We first consider a number of observed properties of the 
local AGN population using our INTEGRAL sample. 

Fig. [3] shows the observed distribution of absorption 
columns for our objects, while Fig. |4] shows the observed 
dependence of the obscured AGN fraction on hard X-ray lu¬ 
minosity. The latter was obtained by counting obscured and 
unobscured sources within specified luminosity bins and di¬ 
viding the first number by the sum of the two. One can 
clearly see a declining trend of the obscured AGN fraction 
with increasing luminosity, which is well known from previ¬ 
ous studies. 

We next calculated the observed hard X-ray LF, (?i(I/obs) 
(number of objects per Mpc^ per log Lobs), of local AGN: 
both in binned and analytic form (see Fig. The analytic 
LF model used throughout this study is a broken power law: 

dNAGN _ A , . 

dlogL “ (L/L.)-!-! + (L/L*)t2 ■ ^ 



Log (cm“^) 

Figure 3. Observed distribution of X-ray absorption columns 
for the INTEGRAL AGN. Unobscured (Ah < 10^2 cm-2), 
lightly obscured (lO^^ + Ah < cm 2) and heavily obscured 
(Ah is 102'* cm“2) objects are shown in blue, magenta and red, 
respectively. 



Figure 4. Observed fraction of obscured (Ah IS 1022 cm~2) 
AGN as a function of observed hard X-ray luminosity for the 
INTEGRAL 7-year survey. 


The binned LF was constructed using the standard l/Knax 
method, whereas the best-fit values (and their uncertainties) 
of the characteristic luminosity, L«, and of the two slopes, 71 
and 72 , of the analytic model (see Tabled w ere found using 
a maximum likelihood estimator (similarly to lSazonov et al.l 
l2007tl : 


£= -2^1n 


0(Z/obs,i) V^ax(-^/obs,i) 

^ 0(-^obs) Vinax{-/^obs)f^log Lohs 


( 2 ) 


© 0000 RAS, MNRAS 000, 000-000 



























Does obscured AGN fraction depend on luminosity? 5 



Log erg s"') 

Figure 5. Observed (in the INTEGRAL 7-year survey) hard X- 
ray luminosity function of local AGN (filled circles) fitted by a 
broken power law (black solid line). The best-fit parameters are 
given in Tabled For comp arison the LF based on the Swift/'QPIT 
survey llAiello^ al.|[2012h is shown by the magenta dashed line. 

where I/obs,i are the observed luminosities of AGN in our 
sample, and V'max(robs) is the volume of the Universe probed 
by the INTEGRAL 7-year survey for a given Lobs, which 
can be calculated from the sky coverage curve (see Fig. [T]). 
The normalization of the analytic model is derived from the 
actual number of objects in the sample. 

Comparing this newl y determined observe d hard X-ray 
LF with our old result (ISazonov et al.l l2007ll based on a 
smaller set (66 vs. 151 objects) of AGN detected with IN¬ 
TEGRAL, we find good agreement between the two, but 
the constraints on the LF parameters have now significantly 
improved. We can also compare the INTEGRAL LF with 
that derived from a still larger (361 objects) sample of 
(mostly) local AGN found in nearly the same energy ban d 
(15-55 keV) in the Swift/HAT survey (lAiello et al.ll2of^ . 
As can be seen in Fig. [S] the two LFs are in good agreement 
with each other. 

Finally, we calculated separately the observed LFs of 
unobscured and obscured AGN (see Fig. El). It can be seen 
that these LFs are different in shape, as is verified by the 
best-fit parameters of the corresponding analytic fits (see 
Table EJ. 


4 INTRINSIC PROPERTIES OF LOCAL AGN 

The observed LF just discussed has not been corrected for 
any effects associated with absorption or scattering of hard 
X-rays emitted by the AGN central source on the way be¬ 
tween the source and the observer, nor for any intrinsic 
anisotropy of the emission generated by the central source in 
AGN. This observed LF is expected to be affected by absorp¬ 
tion bias: an obscured AGN will be inferred to have a lower 
luminosity. Lobs = /obs x dvrL)^ (here, /obs is the measured 
hard X-ray flux and D is the distance to the source), than 



Figure 6. Observed hard X-ray luminosity functions of un¬ 
obscured (Ajj < 10^^ cm”blue filled circles) and obscured 
(Ajj ^ 10^^ cm”red empty squares) AGN, fitted by broken 
power laws (blue dotted and red dashed lines, respectively). The 
best-fit parameters are given in Table [2] 



Figure 7. Torus model. 


its intrinsic (i.e. emitted by the central source) luminosity, 
Lintr, and a source like this can be found in a flux-limited 
hard X-ray survey within a smaller volume of the Uni¬ 
verse than it would be in the absence of X-ray absorption: 
Lmax(Lobs)/Lmax(Lintr) ^ (Lobs/Lintr) ^ . Here, the approx¬ 
imation simbol reflects the fact that AGN obscuration may 
also affect the shape of the measured X-ray spectrum and 
thus the number of photons recorded by a given detector 
with its specific energy response. On the other hand, as dis¬ 
cussed below, unobscured AGN are expected to have higher 
observed luminosities than their intrinsic angular-averaged 
luminosities and can thus be detected within a larger Umax. 
We can correct for both of these biases and obtain an in¬ 
trinsic hard X-ray LF of local AGN. To this end, we use a 
physically motivated obscuration model described below. 
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Table 2. Fits of different hard X-ray luminosity functions by a broken power law 


AGN 

class 

Aagn 

logL* 

71 

72 

10“® Mpc-3 

Num. density 
(logL = 40.5-46.5) 
lO”’* Mpc-3 

Lum. density 
(logL = 40.5-46.5) 
10®® erg s“® Mpc“® 





Observed LF 




All 

151 

43.74 ±0.19 

0.93 ±0.10 

2.33 ±0.15 

1.122 

54 (41 

±82) 

1.57 ±0.20 

Unobscured 

67 

43.98 ± 0.32 

0.83 ±0.18 

2.37 ±0.28 

0.243 

10 (7- 

F25) 

0.45 ± 0.08 

Obscured 

84 

43.65 ±0.21 

0.99 ±0.21 

2.48 ± 0.22 

0.839 

48 (37 

±78) 

1.14 ±0.18 




Intrinsic LF, isotropic emission, 9 = 30° 




All 

151 

43.69 ±0.18 

0.89 ±0.10 

2.36 ±0.16 

1.806 

61 (47 

±96) 

1.97 ±0.25 

Unobscured 

67 

43.87 ±0.31 

0.84 ±0.18 

2.39 ±0.28 

0.201 

7.0 (5.1 

± 16.7) 

0.30 ± 0.05 

Obscured 

84 

43.62 ±0.21 

0.93 ±0.12 

2.45 ±0.21 

1.581 

59 (46 - 

F 101) 

1.63 ±0.26 




Intrinsic LF, isotropic emission, 9 = 45° 




All 

151 

43.71 ±0.19 

0.89 ±0.10 

2.36 ±0.16 

1.494 

52 (41 

±82) 

1.71 ±0.21 

Unobscured 

67 

43.90 ±0.31 

0.84 ±0.18 

2.39 ±0.27 

0.203 

7.5 (5.4 

± 18.2) 

0.32 ± 0.06 

Obscured 

84 

43.62 ±0.21 

0.92 ±0.12 

2.44 ±0.21 

1.698 

59 (46 - 

F 105) 

1.70 ±0.26 




Intrinsic LF, cosine-law emission, 9 = 30° 




All 

151 

43.69 ±0.17 

0.90 ±0.10 

2.43 ±0.17 

1.593 

57 (44 

±87) 

1.79 ±0.21 

Unobscured 

67 

43.70 ± 0.30 

0.86 ±0.17 

2.42 ± 0.28 

0.188 

5.3 (3.9 

± 11.6) 

0.19 ±0.03 

Obscured 

84 

43.68 ±0.21 

0.93 ±0.12 

2.45 ±0.21 

1.565 

66 (51 - 

F 114) 

1.88 ±0.30 


The normalization A is given without an error because this parameter is strongly correlated with the others. 


4.1 Torus model and AGN spectra 

We have a developed a Monte-Carlo code for mod¬ 
elling AGN X-ray spectra modified by reprocessing in 
a toroidal structure of gas. The adopted geometry (see 
Fig. [71) is similar to that used in ot h er existing models, 
e.g. Ilted a. Awak i. fc Terashimd (l2009l l: iMurphv fc YaaooH 
||2009[) : Brightman fc Nandral I 201lli . The key assumptions 
of our model are: 

(i) The geometrical shape is that of a ring torus; 

(ii) The gas is homogeneous, cold, neutral and of normal 
cosmic chemical composition; 

(iii) The X-ray spectrum emitted by the central source 
is a power law with an exponential cutoff, dN/dE oc 

with r = 1.8 and = 200 keV; 

(iv) The central (point-like) source is either isotropic, 
dLinti/dQ = const - hereafter, Model A, or emitting ac¬ 
cording to Lambert’s law, dL^ntr/dQ oc cos a, where a is 
the viewing angle with respect to the axis of the torus - 
hereafter, Model B. 

The introduction of Model B is an important aspect of 
the present study and is motivated by the common belief 
that the hard X-ray emission observed from AGN is pro¬ 
duced by Comptonization of softer emission from an ac¬ 
cretion disk around a SMBH in a hot corona lying above 
the disk. If such a corona has quasi-planar geometry, the 
hard X-ray flux it produces will be collimated along the 
axis of the disk/corona roughly as T oc /r (the exact 
law being dependent on t he photon energy and the op¬ 
tical depth of the corona. IPozdnvakov. Sobol, fc SunvaevI 
1 19831 : ISunvaev fc TitarchukI Il98^ ~l'. where is the cosine 
of the angle between the outgoing direction and the axis 
of the disk/corona. Because the obscuring torus in turn is 
likely coaligned with the accretion disk, the emergent hard 
X-ray radiation will be collimated along the axis of the 
torus. In reality, a signihcant fraction of the coronal emis¬ 
sion is reflected by the underlying accretion disk, but this 
also occurs preferentially along the axis of the disk/torus 


dMagdziarz fc Zdziarskil[l99!Th . Since there is still significant 
uncertainty in the overall physical picture, we introduce a 
simple, energy independent collimation factor dLintr/dfl oc 
cos a to get an idea of how strongly intrinsic collimation of 
hard X-ray emission can affect observed properties of local 
AGN. 

Apart from the two alternatives for the angular depen¬ 
dence of intrinsic emission (Model A or Model B), our model 
has three free parameters: (i) the equatorial column den¬ 
sity, An.eq (the total number of H atoms per cm^ along an 
equatorial line of sight between the central source and the 
observer, (ii) half-opening angle of the torus, 6, and (iii) the 
viewing angle relative to the axis of the torus, a (see Fig. [7]). 

X-ray photons emitted by the central source can scat¬ 
ter multiple times within the torus before they either get 
photoabsorbed in the gas or escape from the system. Our 
radiative t ransfer calculations are based on a method de¬ 
veloped by IChurazov et ^ (l2008lf . The gas in the torus is 
assumed to be neutral, with the relative abundances of all el¬ 
ements as in the Solar photosphere. The following processes 
are included in the simulations: photoelectric absorption, 
Rayleigh and Compton scattering and fluorescence. Pho¬ 
toelectric absor ption is calculated using the data and ap¬ 
proxim ations of IVerner fc Yakovlev! (Il995h and IVerner et al] 
lll99d). For fluore s cence we use the energies and yields from 
iKaastra fc Mew3 (Il993h . Gompton and Rayleigh scattering 
are modelled using di fferential cros s sections provided by 
the GLECS package llKippenl l2004l l of the GEANT code 
dAgostinelli et al]l2003l L Namely, t he Livermore Eval uated 
Photon Data Library (EPDL, see ICullen et al.l ll99Cll l and 
the Klein-Nishina formula for free electrons are used to cal¬ 
culate total cross-sections and the angular distribution of 
scattered photons for each element. 

Pigure[8]shows examples of emergent AGN spectra sim¬ 
ulated using our model of the obscuring torus. As expected, 
for obscured AGN {a > 0), the observed hard X-ray flux can 
be strongly attenuated relative to the emitted flux and for 
high absorption columns (Ah 3 > 10^"* cm“^) the spectrum 
can become reflection-dominated, as the observer will mostly 
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isotropic emission, 6 = 30 °, 



cosine low, 0 = 30°, 

’ ’ H.eq 



Energy (keV) 

Figure 8. Top: Examples of simulated AGN spectra for Model A, 
half-opening torus angle d = 30° and equatorial column density 
Ajj^eq = 10^® cm“^, for various viewing angles. The dashed curve 
shows the intrinsic (angular-averaged) spectrum. The shaded area 
indicates the INTEGRAL/lYilS energy band used for AGN selec¬ 
tion in this work. Bottom: The same, but for Model B. 


see emission reflected from the inner walls of the torus rather 
than emission from the central source transmitted through 
the torus. In this last case, there also appear strong iron 
Ka and K/3 fluorescent lines at 6.4 keV and 7.06 keV. These 
spectral properties and trends for heavily obscured AGN are 
of course well known. 

We further see from Fig. [S] that the spectra observed 
from directions a < 6, corresponding to unobscured AGN, 
also differ from the intrinsic spectrum. Namely, they have 
an excess due to Compton reflection of hard X-rays from the 
torus in the direction of the observer. This hump is located 
approximately within the energy band of 17-60 keV that we 


use for detecting AGN in the INTEGRAL survey. It is ob¬ 
vious that this Compton reflection component should bias 
observed luminosities of unobscured AGN higher in this and 
similar (e.g. Swift/UNT) hard X-ray surveys. Any intrinsic 
collimation of emission along the axis of the obscuring torus 
will make this positive bias even stronger (see the spectrum 
for Model B and a « 0 in the lower panel of Fig. H). This im¬ 
portant aspect is frequently overlooked in AGN population 
studies, even though a reflection component is well known 
to be present in the hard X-ray spectra of unobscured AGN. 


4.2 AGN detection bias 


To quantify biases affecting detection of unobscured and ob¬ 
scured AGN in the INTEGRAL survey, we show in Fig. 
for Model A and Model B, the ratio, R{NH,eq,d,ci) = 
Tobs/Tintr, of the observed to intrinsic luminosity in the 
17-60 keV energy band as a function of An.eq, for a torus 
half-opening angle 6 = 30° and several narrow ranges of the 
viewing angle a. One can see that R is always larger than 
unity, i.e. Lobs > Lintr, for unobscured AGN. For example, 
for Model A and a « 0, 7? increases from 1 to ~ 2 as Nu.eq 
increases to a few 10^“* cm“^ and remains at approximately 
this level thereafter. This trend can be easily understood: 
the amplitude of Compton reflection is expected to be pro¬ 
portional to the torus optical depth, r, in the optically thin 
regime (r ^ 1) and constant in the opposite case (r ^ 1). 
As regards obscured AGN (a > 6), R decreases with increas¬ 
ing Au.eq and increasing a (apart from a local maximum 
at An.eq 3 X 10^^ cm“^ at near-equatorial directions for 
Model B - due to the reflected component), as could be ex¬ 
pected due to the increasing attenuation of the transmitted 
component. The most obvious and important difference of 
Model B with respect to Model A is that the observed hard 
X-ray flux is anisotropic even in the absence of an obscuring 
torus (i.e. for Ah = 0 ) - just due to the initial collimation 
of emission. 

We can proceed further and ask the question: what 
would be the average observed/intrinsic flux ratio for the 
local populations of unobscured and obscured AGN if (i) 
AGN tori were randomly oriented with respect to the ob¬ 
server, which is a natural assumption, and (ii) all the tori 
had the same half-opening angle 6 (this, of course, permits 
the physical size of the torus to vary from one object to an¬ 
other and e.g. to depend on luminosity). To this end, we just 
need to average the dependencies shown in Fig. [9] over the 
viewing angle a for the unobscured and obscured directions: 


77unobsc(AH,eq, ^) — 


iLs e ^(-^H.eq, 0, a)d cos a 


1 — cos 6 


(3) 


Aobsc(AH,eq, d) 


jv:os0 ^('jYH,eq,0,a)dcOSa 

COS 6 


(4) 


The result is shown in Fig.[T0]as a function of An.eq for torus 
half-opening angles 0 — 30°, 45° and 60°. One can see that 
for Model A, Aunobsc (the average observed/intrinsic flux 
ratio for unobscured AGN) reaches a maximum of ~ 1.5-2, 
depending on 6, at An.eq ~ 5 x 10^“^ cm“^, then declines 
to ~ 1.4-1.7 by An.eq ~ 1.5 x 10^® cm“^ and stays at ap¬ 
proximately this level for higher column densities. The aver¬ 
age ratio Aobsc for obscured AGN monotonically decreases 
from 1 to ~ 0.2 as An.eq increases from <C lO^'* cm“^ to 
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isotropic emission, 0 = 30° 




Figure 9. Top: Calculated ratio of observed to intrinsic (angular- 
averaged) luminosity in the 17-60 keV energy band for a torus 
half-opening angle 0 = 30°, for different viewing angles (a), as a 
function of the torus column density for Model A. Bottom: The 
same, but for Model B. 


~ 1.5 X 10^® cm“^ and stays at this level thereafter. In the 
case of a cosine-law emitting source, the situation is qual¬ 
itatively similar, but the contrast between the unobscured 
and obscured directions is more pronounced: it is present 
already at A^h = 0 and increases further, due to Compton 
reflection, with increasing A^h. 

It is obvious from Fig.[T2]that a hard X-ray survey, like 
the ones performed by INTEGRAL and Swift, will hnd un¬ 
obscured AGN more easily than even lightly obscured ob¬ 
jects, let along heavily obscured ones. Our goal now is to 
correct the observed statistical properties of local AGN for 
this obvious bias. 


isotropic emission 



cosine low 



Figure 10. Top: Calculated ratio of observed to intrinsic 
(angular-averaged) luminosity in the 17—60 keV energy band av¬ 
eraged separately over all unobscured and obscured directions 
(q < 0 and Q > 0, respectively) for three values of the torus 
half-opening angle, as a function of the torus column density, for 
Model A. Bottom: The same, but for Model B. 


4.3 Intrinsic distribution of torus column densities 

We can first estimate the intrinsic distribution of the col¬ 
umn densities, AH.eq, of AGN torii. To this end, we need 
to correct the observed Ah distribution (Fig. [S]) for ab¬ 
sorption bias, excluding from the consideration the first. 
Ah < 10^^ cm“^, bin since it pertains to unobscured AGN 
for which our line of sight does not cross the torus. Given 
the fairly small number of obscured AGN, especially of heav¬ 
ily obscured ones, in our sample, we are bound to make 
some simplifying assumptions. For example, we may as¬ 
sume that the intrinsic An.eq distribution does not depend 
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on luminosity. In this case, the intrinsic A^H,eq distribution 
can be estimated simply by dividing the observed one by 
^obL('^H,eq, 0) (and normalising the resulting dependence 
so that its integral over A^h equals unity), with the bias fac¬ 
tor Rohsc = Aobs/Aintr having been discussed in 1)4.21 

In doing this exercise, we assumed that A^u.eq = 
(4/7r)A^H ss 1.27A^h for our obscured AGN. This is because 
we do not know the orientation of our objects apart from 
the fact that some of them are unobscured and hence a < 6, 
while others are obscured and hence a > 6. For our assumed 
torus geometry (see Fig. 0, the line-of-sight column density 
depends on the viewing angle as follows: 


- (SI)'’ 

so that the mean Ah over all obscured directions is 


A] 


Hjobsc — 


r 

Jo 


NH{a)dcos a yr 

n — ~T'^H,eq. 

cos6» 4 


( 6 ) 


Hence the coefficient in the conversion of Ah to An.eq above. 
Note that the Ah values adopted from the literature for 
some of our Compton-thick AGN may already have been as¬ 
cribed the meaning of an equatorial rather than line-of-sight 
column density by the corresponding authors. However, con¬ 
sidering our sample of heavily obscured sources as a whole, 
the information it contains on the absorption columns is 
very heterogeneous, as it is based on various spectral mod¬ 
els used by various authors. Fortunately, a typical expected 
difference between An.eq and Ah for obscured AGN is only 
20 % (see equation ([S)) above) and has negligible impact 
on our results. 

The resulting intrinsic An.eq distribution is presented 
in Fig. 111! It is only weakly dependent on both the as¬ 
sumed half-opening angle 6 of the obscuring torus and the 
assumed emission model (Model A or Model B). This dis¬ 
tribution can be roughly described as log-uniform between 
An.eq = 10^^ and 10^® cm“^, although the upper bound¬ 
ary is, of course, fairly uncertain. A similar result was pre¬ 
viously obt ained using AGN from the Swift /B AT hard X- 
ray survey (iBurlon et al.ll20l'i] : IUeda et al.ll2014lf . Moreover, 
the intrinsic Ah distribution shown in Fig. [TT] is similar 
to the one inferred for optically se lected Seyfert 2 galaxies 
llRisaliti. Maiolino. fc SalvatilllQ^ '). 


4.4 Intrinsic luminosity function 

We now calculate the intrinsic hard X-ray LF of unobscured 
and obscured AGN, (j>{Linti) = dA/d log Lintr. As for the ob¬ 
served LFs discussed in SJH] we use both binned and analytic 
represent at ions. 

For the binned LFs, the procedure is as follows: 

(i) First, based on the observed luminosity Lobs.i and es¬ 
timated torus column density An.eq.i of each source in the 
sample (Table lAlD . we determine its intrinsic hard X-ray 
luminosity as either Lintr.i = Lobs,i/Aunobsc(AH,eq,i, 6*) (for 
unobscured sources) or Lintr.i = Lobs.i/Aobsc(AH,eq,i,6*) (for 
obscured sources), where the ratios Aunobsc and Aobsc are 
calculated as discussed above (from eqs. ([3]) and (|4|, see 
Fig. HOD , assuming some (the same for all objects) torus 
half-opening angle 9 and using Model A or Model B. Here 
again we use the average ratios Aunobsc and Aobsc rather 



cosine low, 0 = 30" 



Log (cm-2) 

Figure 11. Reconstructed intrinsic distribution of column densi¬ 
ties of obscuring tori in local AGN, calculated assuming 9 = 30° 
and either Model A {top) or Model B (bottom). The dotted line 
corresponds to a log-uniform distribution. 


than the viewing-angle dependent i?(AH,eq, 9, a) from which 
they derive for the lack of knowlegde of the orientation of 
our objects. Strictly speaking, this procedure is not fully 
correct, because for given An.eq and 9, a hard X-ray flux 
limited survey will preferentially find objects with smaller 
viewing angles a within the corresponding groups of a < 9 
and a > 9, as is clear from Fig. [Q] However, this may be 
regarded as a next-order correction to the bias considered 
here and does not significantly affect our results, as we verify 
in ^ 

As was said before, the An.eq values for our obscured 
objects are estimated from their measured Ah columns as 
An.eq = 1.27 Ah. However, we cannot determine similarly 
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the torus column densities for our unobscured AGInB There¬ 
fore, we simply assume that A^n.eq = 10^“^ cm“^ for these ob¬ 
jects, since this is approximately the median value of the in¬ 
ferred intrinsic absorption column distribution for obscured 
AGN (see Fig. [TTIl. 

(ii) Second, we calculate for each source the volume of the 
Universe, Vjnax,i(Tobs,i), over which AGN with such observed 
luminosity can be detected in the INTEGRAL survey. Since 
the detection limit for a given hard X-ray instrument (IBIS 
in our case) is actually determined by photon counts, it 
should depend on the observed X-ray spectral shape, which 
for the problem at hand is affected by absorption and reflec¬ 
tion in the torus (see examples of AGN spectra in Fig. El). 

We correct Vmax.i for this effect, but this correction proves 
to be negligible (as is the fc-correction due to cosmological 
redshift). As a result, we obtain essentially the same Vjnax.i 
for our sources as we used in constructing the observed LF 
in 331l- 

(iii) The Hnal step consists of summing up the l/Knax.i 
contributions of the individual sources, i.e. adding the 
1/Umax,i for each AGN of a given class (unobscured or ob¬ 
scured) to the space density of such objects within a lu¬ 
minosity bin containing Lintr.i (rather than Lobs,i) for this 
source. 

To obtain analytic forms of the intrinsic LFs, we use 
the same broken power-law model as for our observed LFs 
but a different likelihood estimator: 

^ 2 \ ^ ^ 0(Lintr,i) J' Unax(Lintr.i, A^H,eq)d log A^H.eq 

J J 0(Lintr)Ljnax(Lintr, Nh ,eq)dl0g Lintrd log An.eq 

Here Lintr.i are the same estimates of the intrinsic lumi¬ 
nosities of our objects as we used before to construct the 
binned intrinsic LFs (i.e. calculated from Lobs.i using the 
actual A’h,! estimates for the obscured AGN and assum¬ 
ing that A^H.eq = 10^“* cm“^ for the unobscured ones), but 
Lmax (Lintr, Nn.eq) is now the volume over which AGN with 
given intrinsic luminosity Lintr and torus column density 
An.eq can be detected in the INTEGRAL survey. To calcu¬ 
late these volumes, we again use the a-averaged quantities 
I?unobsc(^, A’n.eq) (in fitting the intrinsic LF of unobscured 
AGN) and i?obsc(6, A^n.eq) (in fitting the intrinsic LF of ob¬ 
scured AGN). The integrals over dlogAn.eq in equation 0 
are computed from 10^^ to 10^® cm i.e. we assume that 
the intrinsic distribution of torus column densities is log- 
uniform over this range, as suggested by the result of our 
proceeding analysis shown in Fig. 1111 

Figure [T^] shows the resulting intrinsic LFs for unob¬ 
scured and obscured AGN, calculated assuming 6 — 30° for 
Model A and Model B. We see that in the former case, the 
shapes of the intrinsic LFs of unobscured and obscured AGN 
are clearly different from each other, although to a lesser de¬ 
gree that it was for the observed LFs (Fig. [O]) from which 
they derive. However, for Model B the intrinsic LFs of un¬ 
obscured and obscured AGN are not significantly different 
in shape from each other. These conclusions are verified by 

^ In principle, one could try to estimate IVji.eq for unobscured 
AGN from the contribution of the reflection component to the ob¬ 
served spectrum, but that requires high-quality hard X-ray data, 
which is not always available, and is model-dependent. In partic¬ 
ular, the result will depend on the unknown opening angle 0. 


isotropic emission, 9 = 30' 



cosine low, 0 = 30" 



Log (L;„„, erg s"') 

Figure 12. Top: Reconstructed intrinsic hard X-ray LFs of un¬ 
obscured (blue filled circles) and obscured (red empty squares) 
AGN, fitted by a broken power law (the best-fits parameters are 
given in Table Ell (blue dotted line and red dashed line, respec¬ 
tively). Model A is adopted, with 0 = 30°. Bottom: The same, 
but for Model B. 

the best-fit parameters obtained for these LFs (see Table [2]). 
Note that the derived intrinsic LFs (both binned and ana¬ 
lytic ones) are only weakly sensitive to the torus half-opening 
angle 6 that was assumed in constructing them, and nearly 
the same results are obtained for 6 = 30° and 6 = 45°. 
This is due to the weak sensitivity of the i?unobsc and i?obsc 
factors to 6 (see Fig. [iQl). 

The transformation of the observed binned LFs to the 
intrinsic ones can be understood as follows: (i) all unob¬ 
scured AGN making up the LF shift by the same amount, 
logLobs/Lintr ~ 0.1 and 0.3 for Model A and Model B, 
respectively, to the left along the luminosity axis (since 
we have assumed the same equatorial optical depth of the 
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Isotropic emission, 0 = 30° 



cosine low, 0 = 30° 



Log (Li„„, erg s"') 

Figure 13. Top: Reconstructed intrinsic hard X-ray LF of local 
AGN (black filled circles), fitted by a broken power law (black 
solid line, the best-fit parameters are given in Table O. Model A 
is adopted, with 9 = 30°. For comparison, the observed LF is 
also shown (magenta empty squares and dotted line). Bottom: 
The same, but for Model B. 

torus, Nf[,eq = 10^"* cm~^, for all of our unobscured ob¬ 
jects), and (ii) each obscured AGN from the INTEGRAL 
sample moves its own distance to the right-hand side of 
the plot, this shift being small for lightly obscured objects 
(Ah < cm“^) but substantial (up to log Lintr/iobs ~ 1 

for Ah > 10 ^® cm“^) for heavily obscured ones. 

Finally, we can calculate the intrinsic LF of the entire 
local AGN population, by summing up the contributions of 
unobscured and obscured sources. In obtaining the analytic 
fit in this case, we define Vjnax(Lintr, An.eq) as follows: 

Lmax — Vjnax, unobsc(l COS^) -f Vjnax.obsc COS (8) 

where V(nax,unobsc(Lintrj An.eq) aud Vjnax,obsc (Lintr, AH,eq) 


are the corresponding volumes for unobscured and obscured 
AGN. The resulting LF is shown in Fig. [13] and its best-fit 
parameters are presented in Table [2| 

One can see that the total intrinsic LF is not very dif¬ 
ferent from the total observed LF. This means that the two 
effects observed in Fig. 1121 namely the shift of the LF of un¬ 
obscured AGN to lower luminosities and the shift of the LF 
of obscured AGN to higher luminosities almost compensate 
each other, with this conclusion being only weakly sensitive 
to the assumed torus opening angle and angular dependence 
of intrinsic emission. 

4.5 Total AGN space density 

Integration of the total intrinsic and observed LFs over lumi¬ 
nosity suggests that the cumulative hard X-ray luminosity 
density of local AGN may be underestimated by the ob¬ 
served LF by ~ 10-30%, although this increase is statis¬ 
tically insignihcant (see Table m. Specifically, the intrinsic 
luminosity density of AGN with Lintr > 10'^° ® erg s ^ is 
found to be ~ 1.8 x 10®® erg s“® Mpc“® (17-60 keV), with 
the exact value slightly depending on our assumptions (see 
Table O. 

For our assumed intrinsic AGN spectrum {dN/dE oc 
ratio of luminosities in the 2-10 keV 
and 17-60 keV energy bands is about unity. Therefore, the 
luminosity density of AGN with Lintr > 10"‘° ® erg s ® may 
be estimated at ~ 1.8 x 10®® erg s~® Mpc“® also in the stan¬ 
dard X-ray band (2-10 keV). We may compare this value 
with a prediction for 2 = 0 based on a redshift-dependen t in¬ 
trinsic luminosity function derived bv lUeda et alJ (l2014l ) us¬ 
ing a large heterogenous sample of AGN compiled from vari¬ 
ous surveys. Integration of this LF over the luminosity range 
from 10“*® ® to 10"'®'® erg s“® gives 8 x 10®® erg s“® Mpc“® 
(2-10 keV), which is a factor of ~ 2 smaller than the above 
estimate. In reality, the I/intr(17 —60 keV)/I/intr(2 —10 keV) 
ratio may well be 1.5 rather than 1 due to the ex¬ 
pected presence in AGN spectra of a Gompton reflection 
component associated with the accretion disk. In fact, this 
component was already discussed in 1)4.11 as one of the rea¬ 
sons why hard X-ray emission may be intrinsically colli¬ 
mated in AGN and is implicitly taken into account in our 
anisotropic Model B. Taking this spectral component into 
account, we can lower our estimate of the luminosity den¬ 
sity to ~ 1.2 X 10®® erg s~® Mpc~® (2 -10 keV), which is 
still higher than the lUeda et al.l (|2014| ) result by a factor 
of ~ 1.5. The remaining difference may be related to the 
different procedures used in these works to construct the in¬ 
trinsic LFs and to the larger and more complete sample of 
local heavily obscured AGN used in our study. 

4.6 Intrinsic dependence of obscured AGN 
fraction on luminosity 

Similarly to the observed LF, the observed dependence of 
the fraction of obscured AGN on luminosity (Fig.[4| must be 
affected by AGN detection biases. By removing these biases 
one can obtain an intrinsic dependence of the obscured AGN 
fraction on luminosity. To this end, we should simply divide 
the intrinsic LF of obscured AGN by the total intrinsic LF. 
The result is presented in Fig. [14] for 9 — 30°, Model A and 
Model B, and in Fig. [15] for 6 = 45°. 
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Log (L|„,r, erg s'’) 



Log (Li, 


Figure 14. Top: Reconstructed intrinsic fraction of obscured 
AGN as a function of intrinsic hard X-ray luminosity, calculated 
for 9 = 30° and Model A. Bottom: The same but for Model B. 
The dotted lines indicate the fraction of the sky that will be 
screened from the central source by a torus with half-opening 
angle 6» = 20°, 30°, 45° or 60°. 


We see that in the case of isotropic emission the de¬ 
clining trend of obscured AGN fraction with luminosity is 
retained upon removing the absorption bias, although the 
intrinsic obscured fraction at any luminosity is signihcantly 
higher compared to the observed fraction (see Fig. |3]|. We 
can interpret this resnlt in terms of the torus opening an¬ 
gle, i.e. the fraction of the sky that will be shielded from 
the central source by a toroidal structure of gas. To this 
end, we have drawn in Figs. [14] and [15] four horizontal lines 
corresponding to 6 = 20°, 30°, 45° and 60°. We see that 
if the central sources in AGN are isotropic, then the torus 
opening angle must be smaller than 30° in low-luminosity 
objects (Lintr;^10‘*^ ® erg s“^) and increasing to ~ 45°-60° 
in high-luminosity ones (Tintr > 10^^ erg s“^). 

If, however, the emission from the central SMBH is col¬ 
limated as dLintr/dfl oc cos a, then the derived intrinsic de¬ 
pendence of the obscured AGN fraction on luminosity is in 
fact consistent with the opening angle of the torus being 


Figure 15. The same as Fig. 1141 but for 9 = 45°. 

constant with luminosity, namely 9 ~ 30° - see the bottom 
panel in Fig. 1141 

Importantly, these conclusions are almost insensitive to 
the opening angle (within the range 6 ~ 20-60°) of the torus 
that we actually assumed in deriving the intrinsic luminosity 
dependences of the obscured AGN fraction (compare the 
results for 6 = 45° and 6 — 30° in Figs. [14] and [T5]l . This 
is again due to the fact that the /Junobsc and i?obsc ratios, 
which characterise observational biases for our AGN sample 
and which we have corrected for, are, for a given emission law 
(Model A or Model B), primarily determined by the torus 
column density (see Fig. [IS- Hence, the derived intrinsic 
luminosity dependences of the obscured AGN fraction are 
quite robust but depend on the actual degree of collimation 
of the AGN central source. 


5 DIRECT CONVOLUTION MODEL 

We have demonstrated that positive bias with respect to un¬ 
obscured AGN and negative bias with respect to obscured 
AGN in flux-limited hard X-ray surveys together strongly af¬ 
fect the observed dependence of the obscured AGN fraction 
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on luminosity. Our preceeding analysis consisted of conver- 
ing the observed LFs of unobscured and obscured AGN to 
the intrinsic LFs of these populations. In solving this ’inverse 
problem’, we used a number of simplifications that we noted 
were unlikely to have significant impact on our results. In 
particular, we used the viewing angle-averaged conversion 
factors I^unobsc(Lintr , A^H,eq) and I^obsc (Lintr , A^H,eq) rather 

than the I?unobsc(Lintr, A^H,eq, rr) and I?obsc(Lintr, A^H,eq, rr) 

ratios from which they derive. We also assumed a fixed torus 
column density, A^n.eq = 10^^ cm“^, for all of our unob¬ 
scured AGN. To verify that these assumptions were reason¬ 
able, we now perform a ’direct convolution’ test, as described 
below: 

(i) Assume that the tori in local AGN have the same half¬ 
opening angle 6. 

(ii) Assume that local AGN are oriented randomly with 
respect to us. 

(iii) Assume that the AGN central source is isotropic 
or, alternatively, emitting according to Lambert’s law 
(dl/intr/dfl oc cos a). 

(iv) Assume that the intrinsic distribution of AGN torus 
column densities A^n.eq does not depend on luminosity and 
is log-uniform between A^H.min and A^n.max. Such a distribu¬ 
tion, with Nn.min 10^^ cm“^ and A’H.max ~ 10^® cm“^, 
approximately matches the real NH.eq distribution we have 
inferred using the INTEGRAL sample (see Fig. [TT1) . 

(v) Adopt the intrinsic AGN LF as derived in our pre¬ 
ceeding analysis for given 9 and emission law (see Table [S]). 

(vi) Use the above set of assumptions specifying the in¬ 
trinsic properties of the local AGN population to simulate, 
using our torus obscuration model, AGN properties as would 
be observed in the INTEGRAL survey. 

The main difference with respect to the inverse prob¬ 
lem is that the Lobs/Lintr ratio now explicitly depends on 
the viewing angle, which is randomly drawn for each simu¬ 
lated source, and on A^n.eq, which is drawn from the assumed 
log-uniform distribution for each simulated source, both for 
obscured AGN and for unobscured ones. 

Figure [16] shows the simulated luminosity dependences 
of the observed obscured AGN fraction for 9 = 30° and 
either isotropic or cosine-law emission; Fig. [17] shows the 
corresponding results for 9 = 45°. As our baseline Nu.eq 
distribution we use A^H.min = 10^^ cm“^ and A^n.max = 
10^® cm“^ (solid lines), but we also show results obtained 
for Ah. min = 10^^ cm“^ and Ah. max = 10^® cm“^ and for 
Ah. min = 10^^ cm“^ and Ah. max = 10^® cm“^. The results 
of simulations are compared with the luminosity dependence 
actually observed with INTEGRAL. 

We see that the luminosity dependence of the observed 
obscured AGN fraction predicted for the case of isotropic 
emission and 9 = 30° is inconsistent with the INTEGRAL 
data (x^ = 32.6 per 8 data points between Lobs = 10^^ ® 
and 10'^®'® erg s“^). The isotropic model with a larger torus 
opening angle, 9 = 45°, provides a better match but the fit 
is nevertheless poor (x^ = 15.2 per 8 data points). Among 
the four presented cases, the best agreement between simu¬ 
lations and observations is achieved in the case of cosine-law 
emission and 9 = 30° (x^ = 9-0 per 8 data points). These 
results confirm our previously reached conclusion that un¬ 
less hard X-ray emission in AGN is intrinsically collimated. 




Log erg s"’) 

Figure 16. Simulated dependence of the observed fraction of ob¬ 
scured AGN on observed hard X-ray luminosity for 9 = 30° and 
either isotropic {top) or cosine-law (bottom) emission, for various 
ranges of An.eq (solid and dotted lines). The intrinsic obscured 
fraction is indicated by the dashed line. For comparison, the cor¬ 
responding dependence observed with INTEGRAL is reproduced 
from Fig.(data points with error bars). 


there must be an intrinsic declining trend of the torus open¬ 
ing angle with increasing AGN luminosity. 

Nevertheless, Figs. [T6l and 1171 clearly demonstrate that 
biases associated with detection of AGN in flux-limited hard 
X-ray surveys inevitably lead to the observed fraction of 
obscured AGN being dependent on luminosity even if this 
quantity has no intrinsic luminosity dependence. Specifically, 
in such a case, the observed obscured fraction approaches 
a constant value in the limit of L L*, and another, 
lower limiting value at L ^ L,. It is easy to show that, 
if the intrinsic ratio of obscured and unobscured AGN is 
(Aobsc/Aunc>bsc)intr, then their observed ratio will be: 
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where (-Robac) and (Runobsc) are the appropriately ensemble- 
averaged bias factors (= Robs/Rintr) for unobscured and ob¬ 
scured AGN, respectively, and 7 is the (effective) slope of 
the luminosity function. For example, in the case of cosine- 
law emission and 8 = 30°, the corresponding bias fac¬ 
tors averaged over the log-uniform A^n.eq distribution (see 
Fig. [ini) are (Runobsc) ~ 2 and (Robac) 0.7, whereas 
(Aobac/A'unobac)intr = 6.46. Therefore, for our inferred intrin¬ 
sic AGN LF, with 7 ~ 0.9 and 7 ~ 2.4 in the low- and high- 
luminosity ends, respectively, and L* ~ we may ex¬ 

pect (Aobac/Nunobac)oba 2.5 and (Aobac/Nunobac)oba ^ 0.5 
at L -C 10'^® ’’ and L S> 10“*®'^ erg s“^, respectively. This cor¬ 
responds to the obscured AGN fractions Aobac/(Nunobac + 
Aobac) ~ 0.71 and ~ 0.33, respectively, which is approxi¬ 
mately what we see in Fig. [16] for the results of simulations 
with the cosine-law emission law and 6 — 30°. 

Figure [18] shows how our direct convolution model pre¬ 
dicts the observed LFs of unobscured and obscured AGN in 
the case of cosine-law emission and 6 — 30°. We see that for 
unobscured AGN, the simulated observed LF fits the INTE- 


cosine low, 0 = 30" 



cosine low, 0 = 30" 



Log erg s"') 

Figure 18. Top: Simulated observed LF of unobscured AGN 
(dashed line) for 0 = 30°, cosine-law emission, AH.min = 
10^^ cm“^ and Ajj.max = 10^® cm“^, in comparison with the 
LF of unobscured AGN observed by INTEGRAL (data points 
with error bars and their fit by a broken power-law model — dotted 
line). The solid line shows the assumed intrinsic LF of unobscured 
AGN (equal to the total AGN LF multiplied by (1 — cos 0)), which 
was previously (in il4.4ll derived from the INTEGRAL data using 
an inverse approach. Bottom: The same, but for obscured AGN. 
In this case the intrinsic LF of obscured AGN is equal to the total 
intrinsic LF multiplied by cos 0. 


ORAL data well. In the case of obscured AGN, the match 
between the simulated and actually observed LFs is good 
below Lobs ~ erg s“^, but a significant deviation is 

evident at higher luminosities. This probably reflects the in¬ 
trinsic differences between the inverse and direct approaches 
to the considered problem, discussed at the beginning of this 
section. Nevertheless, this difference does not significantly 
affect the conclusions of this study. 

Note that in reality the situation may be more compli¬ 
cated. For example, the intrinsic Nu.eq distribution may de¬ 
pend on luminosity. Also, the obscuring gas may be clumpy. 
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so that the A^h distribution observed for obscured AGN may 
represent not only the distribution of torus column densities 
over the AGN population but also the distribution of line-of- 
sight columns over different observing directions for a given 
AGN torus. It would be interesting to study this and other 
possibilities in future work, when significantly larger samples 
of hard X-ray selected AGN become available. 


6 DISCUSSION AND SUMMARY 

We utilised a sample of about 150 local («<0.2) hard X-ray 
selected AGN, with reliable information on X-ray absorption 
columns, to find out how strongly the observed declining 
trend of the obscured AGN fraction with increasing lumi¬ 
nosity may be affected by selection effects. Using a torus- 
obscuration model and a state-of-the-art radiative transfer 
code, we demonstrated that there must exist not only a neg¬ 
ative bias, due to absorption in the torus, in finding obscured 
AGN in hard X-ray flux limited surveys, but also a positive 
bias in detecting unobscured AGN - due to reflection by the 
torus of part of the radiation emitted by the central source 
towards the observer. We further pointed out that these two 
biases may in fact be even stronger if one takes into account 
plausible intrinsic collimation of hard X-ray emission along 
the axis of the obscuring torus, which can arise both in the 
hot corona where the hard X-ray emission presumably orig¬ 
inates and as a result of reflection of part of this radiation 
by the underlying, optically thick accretion disk. 

We demonstrated that for an AGN luminosity function 
that steepens at high luminosities, which is indeed the case, 
these observational biases should inevitably lead to the ob¬ 
served fraction of obscured AGN being smaller in the high- 
luminosity end of the LF than in the low-luminosity end 
even if the obscured AGN fraction has no intrinsic luminos¬ 
ity dependence. Moreover, even in the low-luminosity part 
of the LF, the observed obscured fraction will be lower than 
its intrinsic value. 

We explored two possibilities for the central hard X- 
ray source in AGN: (i) isotropic emission and (ii) emission 
collimated according to Lambert’s law, dL/dQ, oc cos a. In 
the former case, the intrinsic (i.e. corrected for the biases 
discussed above) obscured AGN fraction reconstructed from 
our INTEGRAL sample still shows a declining trend with lu¬ 
minosity, although the inferred intrinsic obscured fraction is 
larger than the observed one at any luminosity. Namely, the 
obscured fraction is larger than 85% at erg s“^ 

(17-60 keV), and decreases to <60% at L > erg s“^. In 
terms of the half-opening angle 6 of an obscuring torus, this 
implies that ^;^30° in lower-luminosity AGN, and 9 ^ 45° 
in higher-luminosity ones. If, however, the emission from the 
central SMBH is collimated as dL/dO, oc cosq, then the de¬ 
rived intrinsic dependence of the obscured AGN fraction is 
consistent with the opening angle of the torus being constant 
with luminosity, namely 6 30°. 

At the moment, we regard both possiblities - intrin¬ 
sic obscuring AGN fraction declining with luminosity or be¬ 
ing constant - as feasible, as they depend on the presently 
poorly understood angular emission diagram of the cen¬ 
tral source in AGN. We note however that a luminosity- 
independent obscured AGN fraction might be consistent 
with findings of some studies based on non-X-ray selected 


AGN samples (see [Lawrence fc ElvT3l2010l for a discussion). 
We also note that the intrinsic ratio of obscured to unob¬ 
scured AGN that follows from our study, which changes with 
luminosity from > 6 : 1 to ~ 1 : 1 in the case of isotropic 
emission and is ~ 6 : 1 in the case of cosine-law emission, 
is not very different fro m the ~ 4 : 1 ratio inferr ed for op¬ 
tically selected AGN by iMaiolino fc Riekd (ll995l L A more 
careful comparison of these and other existing estimates of 
the ratio of obscured and unobscured AGN in future work 
may help us get insight into the geometrical and physical 
properties of obscuration in AGN, which may be different in 
X-ray, optical, infrared and radio bands. 

The intrinsic dependence of the obscured AGN frac¬ 
tion on luminosity derived here can find application in 
modelling the cosmic X-ray background (CXB). Impor¬ 
tantly, the inferred obscured fractions are somewhat larger 
(even without allowance for possible intrinsic collimation 
of X-ray emission in AGN) than t hose adopted in some 
popular CXB synthesis models (e .g. TLeister&^Jr^ l2005l : 
iGilli. Comastri. fc Hasingeil 120071 : lUeda et al.l 2014) - for 
example compare our Figs. fTT and [T^ with Fig. 13 in 


Gilli. Comastri. fc Hasingerl 12007 and Fig. 5 in lUeda et al.l 
20141 (note, however, that these plots use the 2-10 keV en¬ 


ergy band while we use 17-60 keV). 

As a byproduct, we reconstructed the intrinsic hard X- 
ray luminosity function of local AGN and estimated the total 
number density and luminosity density of AGN with L > 
erg s“^ (17-60 keV and 2-10 keV), which may be 
used as reference 2 = 0 values in the study of cosmic AGN 
evolution and in modelling the cosmic X-ray background. 

The constraints on the intrinsic dependence of the ob¬ 
scured AGN fraction on luminosity obtained in this work can 
be improved in the near future using larger samples of hard 
X-ray selected AGN from INTEGRAL, Swift and NuSTAR 
surveys Isee lLansburv et al.ll2015l for a new constraint on the 
abundance of heavily obscured AGN from NuSTAR data). 
Note however that it will be practially impossible to improve 
the current, fairly uncertain estimate of the obscured AGN 
fraction at the highest luminosities (> 10'^® erg s“^) in the 
local Universe, since the INTEGRAL and Swift all-sky sur¬ 
veys are sensitive enough to detect all such objects in the 
local (2<0.2) Universe and have found just a few of them 
because of the very low space density thereof. 

We finally note that AGN selection effects similar to 
those discussed for hard X-ray surveys in this work should 
also affect obscured AGN fractions inferred from samples 
selected in X-rays (at energies below 10 keV). This should 
be studied in future work. 
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APPENDIX A: AGN GATALOGUE 

Table ED presents the sample of non-blazar AGN at 
|6| > 5° used in this work. It is bas ed on the cat¬ 
alogue of sources iKrivonos et al.l l2010bl l detected dur- 
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ing the INTEGRAL/IBIS 7-year all-sky survey. After 
publication of this catalogue, three previously unidenti¬ 
fied sources, IGR J13466-1-1921, IGR J14488—4009 and 
IGR J17036-1-3734 have been proved to be AGN and hence 
added to our sample. For 31 nearby (closer than 40 Mpc) 
objects we adopt distance estimates from the Extragalactic 
Distance Database (EDD (0, whereas the distances of the re¬ 
maining objects are estimated from their redshifts, which are 
adopted from NED. The quoted hard X-ray luminosities are 
observed ones, calculated from the adopte d distances and 
measu red hard X-ray (17-60 keV) fluxes (iKrivonos et al.l 

l2010bll . 

For this study, the most important AGN property is 
X-ray absorption column density. The corresponding infor¬ 
mation h as been updated with respect to our previous pub¬ 
licati ons (ISazonov fc RevnivtsevI [20041 : l^zonov et al.ll200'4 . 
l2012|j whenever necessary and possible. The last column of 
Table [AT] provides relevant references. 
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Table Al. INTEGRAL 7-year sample of non-blazar AGN at |fe| > 5°. 


Object z D Ref. logLi 7 _eo keV Ref. 




Mpc 


erg s ^ 

10^^ cm“^ 


IGR J00040+7020 

0.0960 

442.4 


44.31 

3 

2 

SWIFT J0025.8-I-6818 

0.0120 

52.0 


42.51 

> 1000 

3 

MRK 348 

0.0150 

65.2 


43.69 

30 

4 

ESO 297-G018 

0.0252 

110.4 


43.86 

50 

5 

IGR J01528-0326 

0.0172 

74.9 


43.06 

14 

6 

NGC 788 

0.0136 

59.0 


43.33 

40 

4 

MRK 1018 

0.0424 

188.1 


43.69 

< 1 

7,8 

IGR J02086-1742 

0.1290 

607.6 


44.77 

< 1 

9 

LEDA 138501 

0.0492 

219.4 


44.18 

< 1 

10 

MRK 590 

0.0264 

115.7 


43.26 

< 1 

11 

SWIFT J0216.3-I-5128 

0.4220 

2322.8 


46.15 

3 

7 

MRK 1040 

0.0167 

72.7 


43.40 

< 1 

12 

IGR J02343-I-3229 

0.0162 

70.5 


43.26 

2 

13 

NGG 985 

0.0431 

191.3 


43.91 

< 1 

14 

NGG 1052 

0.0050 

19.4 

1 

41.89 

20 

15 

NGG 1068 

0.0038 

12.3 

1 

41.58 

> 1000 

16 

IGR J02524-0829 

0.0168 

73.1 


42.97 

12 

9 

NGG 1142 

0.0288 

126.5 


43.99 

50 

4 

NGG 1194 

0.0136 

59.0 


42.82 

~ 100? 

17 

IGR J03249-I-4041 

0.0476 

212.0 


43.84 

3 

18 

IGR J03334-I-3718 

0.0550 

246.3 


44.17 

< 1 

7,18 

NGG 1365 

0.0055 

18.0 

1 

42.16 

~ 50 

19 

ESO 548-G081 

0.0145 

63.0 


43.21 

< 1 

20 

3C 105 

0.0890 

408.2 


44.69 

30 

7,8 

3G 111 

0.0485 

216.1 


44.53 

< 1 

21 

IRAS 04210-1-0400 

0.0450 

200.0 


43.98 

30 

7,18 

3C 120 

0.0330 

145.4 


44.14 

< 1 

11 

UGG 03142 

0.0217 

94.8 


43.59 

3 

7 

GGGG 420-015 

0.0294 

129.2 


43.51 

> 1000 

3 

ESO 033-G002 

0.0181 

78.9 


43.05 

1 

22 

LEDA 075258 

0.0160 

69.6 


42.75 

< 1 

7,18 

XSS J05054-2348 

0.0350 

154.4 


44.13 

6 

23 

IRAS 05078-1-1626 

0.0179 

78.0 


43.66 

< 1 

24 

ARK 120 

0.0327 

144.0 


44.02 

< 1 

25 

ESO 362-G018 

0.0124 

53.8 


43.10 

< 1 

26 

PIG A 

0.0351 

154.9 


44.01 

< 1 

11 

NGG 2110 

0.0078 

29.0 

1 

43.11 

14 

27 

MCG 8-11-11 

0.0204 

89.0 


43.92 

< 1 

11 

MRK 3 

0.0135 

58.6 


43.45 

100 

28 

PMN J0623-6436 

0.1289 

607.1 


44.73 

< 1 

29 

IGR J06239-6052 

0.0405 

179.4 


43.80 

20 

30 

IGR J06415-I-3251 

0.0172 

74.9 


43.39 

16 

31 

MRK 6 

0.0188 

81.9 


43.52 

~ 5 

32 

IGR J07563-4137 

0.0210 

91.7 


43.08 

< 1 

33 

ESO 209-G012 

0.0405 

179.4 


43.80 

< 1 

34 

IGR J08557-I-6420 

0.0370 

163.5 


43.64 

20 

7,18 

IRAS 09149-6206 

0.0573 

257.0 


44.19 

< 1 

10 

IGR J09253-I-6929 

0.0390 

172.6 


43.66 

> 100? 

35 

NGG 2992 

0.0077 

29.0 

1 

42.75 

1 

36 

MCG -5-23-16 

0.0085 

36.8 


43.27 

2 

11 

IGR J09522-6231 

0.2520 

1276.8 


45.37 

6 

37 

NGC 3081 

0.0080 

28.6 

1 

42.64 

100 

38 

ESO 263-G013 

0.0333 

146.7 


43.72 

30 

39 

NGC 3227 

0.0039 

26.4 

1 

42.84 

< 1 

4 

NGC 3281 

0.0107 

46.3 


43.08 

200 

40 

IGR J10386-4947 

0.0600 

269.6 


44.05 

1 

4 

IGR J10404-4625 

0.0239 

104.6 


43.48 

3 

4 

NGC 3516 

0.0088 

38.0 

1 

42.87 

< 1 

11 

NGC 3783 

0.0097 

25.1 

1 

42.95 

< 1 

4 

IGR J11459-6955 

0.2440 

1230.9 


45.31 

< 1 

7,41 

IGR J12009-I-0648 

0.0360 

159.0 


43.67 

11 

31 
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Table Al. (continued) 


Object z D Ref. logLi 7 _ 60 keV Ref- 




Mpc 


erg s ^ 

10^^ cm“2 


IGR J12026-5349 

0.0280 

122.9 


43.75 

2 

33 

NGC 4051 

0.0024 

17.1 

1 

41.97 

< 1 

11 

NGC 4138 

0.0030 

13.8 

1 

41.60 

8 

42 

NGC 4151 

0.0033 

11.2 

1 

42.66 

8 

11 

IGR J12107+3822 

0.0229 

100.1 


43.14 

3 

7 

NGC 4235 

0.0080 

31.5 

1 

42.06 

< 1 

7 

NGC 4253 

0.0129 

56.0 


42.74 

< 1 

43 

NGC 4258 

0.0015 

7.6 

1 

40.95 

7 

44 

MRK 50 

0.0234 

102.4 


43.21 

< 1 

24 

NGC 4388 

0.0084 

16.8 

1 

42.78 

40 

36,45 

NGC 4395 

0.0011 

4.7 

1 

40.55 

2 

46 

NGC 4507 

0.0118 

51.2 


43.56 

~ 70 

47 

ESO 506-G027 

0.0250 

109.5 


43.76 

100 

48 

XSS J12389-1614 

0.0367 

162.1 


43.98 

2 

33 

NGC 4593 

0.0090 

37.3 

1 

42.96 

< 1 

4 

WKK 1263 

0.0244 

106.8 


43.39 

< 1 

31 

NGC 4939 

0.0104 

34.7 

1 

42.36 

> 1000? 

49 

NGC 4945 

0.0019 

3.4 

1 

41.41 

400 

50 

ESO 323-G077 

0.0150 

65.2 


43.06 

30 

4 

IGR J13091+1137 

0.0251 

109.9 


43.64 

60 

39 

IGR J13109-5552 

0.1040 

481.8 


44.71 

< 1 

7 

IGR J13149+4422 

0.0366 

161.7 


43.68 

5 

13 

MCG -03-34-064 

0.0165 

71.8 


43.19 

~ 50 

36 

CEN A 

0.0018 

3.6 

1 

41.99 

11 

27 

ESO 383-G018 

0.0124 

53.8 


42.74 

20 

7,51 

MCG -6-30-15 

0.0077 

25.5 

1 

42.48 

< 1 

12 

NGC 5252 

0.0230 

100.6 


43.91 

5 

36 

MRK 268 

0.0399 

176.7 


43.81 

30 

7,18 

IGR J13466-I-1921 

0.0850 

388.7 


44.62 

< 1 

52 

IC 4329A 

0.0160 

69.6 


44.02 

< 1 

12 

LEDA 49418 

0.0509 

227.2 


43.75 

2 

53 

NGC 5506 

0.0062 

21.7 

1 

42.92 

3 

11 

IGR J14175-4641 

0.0766 

348.3 


44.21 

> 100? 

35 

NGC 5548 

0.0172 

74.9 


43.27 

< 1 

11 

ESO 511-G030 

0.0224 

97.9 


43.45 

< 1 

11 

NGC 5643 

0.0040 

11.8 

1 

41.23 

> 1000 

3 

NGC 5728 

0.0094 

24.8 

1 

42.51 

200 

3 

IGR J14488-4009 

0.1230 

577.1 


44.69 

6 

54 

IGR J14552-5133 

0.0160 

69.6 


42.87 

< 1 

10 

IGR J14561-3738 

0.0246 

107.7 


43.20 

~ 100 

37 

IC 4518A 

0.0157 

68.3 


43.07 

10 

7,18 

MRK 841 

0.0364 

160.8 


43.98 

< 1 

55 

NGC 5995 

0.0252 

110.4 


43.66 

< 1 

11 

IGR J15539-6142 

0.0149 

64.8 


42.60 

20 

10 

ESO 389-G002 

0.0194 

84.6 


42.93 

6 

18 

WKK 6092 

0.0156 

67.8 


42.96 

< 1 

24 

IGR J16185-5928 

0.0350 

154.4 


43.55 

10 

56 

ESO 137-G034 

0.0092 

33.0 

1 

42.26 

> 1000 

39 

IGR J16385-2057 

0.0264 

115.7 


43.27 

< 1 

13 

IGR J16482-3036 

0.0313 

137.7 


43.85 

< 1 

4 

NGC 6221 

0.0050 

15.6 

1 

41.66 

1 

57 

NGC 6240 

0.0245 

107.3 


43.76 

250 

3 

IGR J16558-5203 

0.0540 

241.6 


44.22 

< 1 

4 

IGR J17009-I-3559 

0.1130 

526.7 


44.74 

30 

7,18 

IGR J17036-I-3734 

0.0650 

293.1 


44.36 

< 1 

18 

NGC 6300 

0.0037 

13.1 

1 

42.02 

25 

11 

IGR J17418-1212 

0.0372 

164.4 


43.86 

< 1 

4 

H 1821+643 

0.2970 

1541.0 


45.58 

< 1 

11 

IC 4709 

0.0169 

73.6 


43.36 

12 

23 

IGR J18249-3243 

0.3550 

1895.4 


45.61 

< 1 

58 

ESO 103-G035 

0.0133 

57.7 


43.43 

30 

11 
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Table Al. (continued) 


Object 

z 

D 

Mpc 

Ref. 

logLi7_60 keV 
erg s~^ 

Ah 

10^^ cm“2 

Ref. 

3C 390.3 

0.0561 

251.4 


44.59 

< 1 

4 

ESO 140-G043 

0.0142 

61.7 


43.14 

2 

59 

ESO 025-G002 

0.0289 

126.9 


43.50 

< 1 

7,18 

IGR J18559+1535 

0.0838 

382.9 


44.54 

< 1 

4 

EXSS 1849.4-7831 

0.0420 

186.3 


44.02 

< 1 

18 

IGR J19077-3925 

0.0760 

345.4 


44.21 

< 1 

7,18 

IGR J19194-2956 

0.1668 

804.6 


45.03 

< 1 

7,18 

ESO 141-G055 

0.0371 

164.0 


44.06 

< 1 

11 

SWIFT J1930.5-I-3414 

0.0633 

285.1 


44.10 

30 

6 

IH 1934-063 

0.0106 

45.9 


42.56 

< 1 

60 

IGR J19405-3016 

0.0520 

232.3 


43.98 

< 1 

61 

NGC 6814 

0.0052 

22.0 

1 

42.42 

< 1 

12 

XSS J19459+4508 

0.0539 

241.1 


44.03 

11 

33 

CYG A 

0.0561 

251.4 


44.67 

20 

62 

ESO 399-IG020 

0.0250 

109.5 


43.25 

< 1 

7,56 

IGR J20286+2544 

0.0142 

61.7 


43.21 

50 

2,4 

4G +74.26 

0.1040 

481.8 


44.99 

< 1 

11 

MRK 509 

0.0344 

151.7 


44.21 

< 1 

11 

RX J2044.0+2833 

0.0500 

223.1 


44.08 

< 1 

7,18 

S5 2116+81 

0.0860 

393.6 


44.72 

< 1 

63 

IGR J21196+3333 

0.0510 

227.7 


43.95 

< 1 

7,18 

NGC 7172 

0.0087 

31.9 

1 

42.88 

13 

11 

IGR J22292+6646 

0.1120 

521.7 


44.54 

< 1 

58 

NGC 7314 

0.0048 

15.9 

1 

41.92 

1 

64 

MRK 915 

0.0241 

105.5 


43.45 

3 

7,18 

MR 2251-178 

0.0640 

288.4 


44.72 

< 1 

65 

NGC 7465 

0.0066 

26.5 

1 

42.32 

~ 10 

66 

NGC 7469 

0.0163 

70.9 


43.44 

< 1 

4 

MRK 926 

0.0469 

208.8 


44.28 

< 1 

4 
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